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We use the framework set up recently to compute non-perturbatively inelastic scattering from 
quantum impurities [G. Zarand et al, Phys. Rev. Lett. 93, 107204 (2004)] to study the the 
energy dependence of the single particle S-matrix and the inelastic scattering cross section for a 
number of quantum impurity models. We study the case of the spin S — 1/2 two-channel Kondo 
model, the Anderson model, and the usual S = 1/2 single-channel Kondo model. We discuss the 
difference between non-Fermi liquid and Fermi liquid models and study how a cross-over between 
the non- Fermi liquid and Fermi liquid regimes appears in case of channel anisotropy for the 5=1/2 
two-channel Kondo model. We show that for the most elementary non-Fermi liquid system, the 
two-channel Kondo model, half of the scattering remains inelastic even at the Fermi energy. Details 
of the derivation of the reduction formulas and a simple path integral approach to connect the 
T-matrix to local correlation functions are also presented. 

PACS numbers: 75.20.Hr,74.70.-b 



I. INTRODUCTION 

Quantum interference effects play a major role in 
mesoscopic systems: they lead to phenomena such 
as Aharonov-Bohm interference^, weak localization ef- 
fects 2 -^, universal conductance fluctuations^, or meso- 
scopic local density of states fluctuations 5 -. All these in- 
teresting phenomena rely on the phase coherence of the 
conduction electrons. This phase coherence is, however, 
destroyed through inelastic scattering processes, where 
some excitation is created in the environment, and which 
thus lead to a loss of quantum interference after a char- 
acteristic time, Tip. This characteristic time is called the 
dephasing time or sometimes as the inelastic scattering 
time. The excitations created in course of an inelastic 
scattering process may be phonons, magnons, electro- 
magnetic radiation, or simply electron-hole excitations, 
where in the latter case the 'environment' is provided by 
the conduction electrons themselves. 

A few years ago Mohanty and Webb measured the de- 
phasing time t v (T) very carefully down to very low tem- 
peratures through weak localization experiments, and re- 
ported a surprising saturation of it at the lowest tempera- 
tures£i£ These experiments gave rise to many theoretical 
speculations: intrinsic dephasing due to electron-electron 
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scattering from two-level systems 
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interaction 

and even coupling to zero point fluctuations have been 
proposed to explain the observed saturation, and induced 
rather harsh discussions. Finally, it has been shown re- 
cently that an apparent saturation can also appear due 
to inelastic scattering from magnetic impurities i&i 3 - 

Triggered by these results of Mohanty and Webb, a 
number of experimental groups also revisited the problem 
of inelastic scattering and dephasing in quantum wires 



and disordered metals: A series of experiments have been 
performed to study the non-equilibrium relaxation of the 
energy distribution function in short quantum wires of 
various compositions. 14 Depending on the material, these 
energy relaxation experiments could be well explained in 
terms of the orthodox theory of electron-electron interac- 
tion in one-dimensional wires ^ and/or inelastic scatter- 
ing mediated by magnetic impurities i^ i 16 i 17 i 18 i 19 i 20 Par- 
allel to, and partially triggered by these experiments, a 
systematic study of the inelastic scattering from mag- 
netic impurities has also been carried out recently, where 
inelastic scattering at energies down to well below the 
Kondo scale has also been studied i 21 ' 22 > 23 Describing in- 
elastic scattering from magnetic impurities around and 
below the Kondo scale has been a theoretical challenge, 
since this regime can be reached only through non- 
perturbative methods. This goal has been finally reached 
in Refs. [H and [24]: In Ref. [TH a theory of inelastic scat- 
tering has been developed at T = temperature, while 
the authors of Ref. [2J showed that the finite temper- 
ature version of the formula introduced in Ref. |l3| de- 
scribes the dephasing rate that appears in the expression 
of weak localization in the limit of small concentrations 
too. Except for very low temperatures, where a small 
residual inelastic scattering is observed ) 22 i 23 these calcu- 
lations are in very good agreement with the experiments 
that clearly show that magnetic impurities in concentra- 
tion as small as lppm can induce substantial inelastic 
scattering j 13 i 24 ' 25 The source of the small residual inelas- 
tic scattering is unclear: It may be due to some mis- 
positioned magnetic impurities with anomalously small 
Kondo temperature or structural defects created by the 
ion implantation, but an intrinsic effect cannot be ex- 
cluded either, although the residual dephasing seems to 
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be proportional to the impurity concentration. Further- 
more, we have to emphasize, that other experiments on 
very dirty metals probably cannot be explained in terms 
of magnetic scattering, and possibly other mechanisms 
are needed to account for the dephasing observed at very 
low temperatures in these systems^ 

The purpose of the present paper is to demonstrate, 
how the rather general theory of Ref. [IH can be applied 
to various quantum impurity problems. In our previ- 
ous work we presented results only for the single channel 
5 = 1/2 Kondo model, while in the present paper we 
extend our study to different quantum impurity models 
(two channel 5=1/2 Kondo model, Anderson model) 
as well, and we also discuss some analytical expressions 
for various scattering rates in the single channel Kondo 
model. In addition, we present many details of derivation 
of the formalism shortly discussed in Ref. [l3|. 

Ref. [TH formulates the problem of inelastic scattering 
in terms of the many-body 5-matrix defined through the 
overlap of incoming and outgoing scattering states: 



out (/I*) in = ln(f\S\i)i 



(1) 



The incoming and outgoing scattering states, |i)i n and 
|/)out, are asymptotically free, and they may contain 
many excitations, i.e. they are true many-body states. 
In the interaction representation 5 is given by the well- 
known expression 



5 = Texp 



H int (t)dt 



(2) 



with T the usual time-ordering operator, and H{ nt the 
interaction part that induces scattering. 

The many-body T-matrix is defined as the 'scattering 
part' of the 5-matrix, 



S = I + iT, 



(3) 



where / denotes the identity operator. Energy conserva- 
tion implies that 



l (f\T\i) ia = 2*8(E f -E i ) (f\T\i) 



(4) 



with the (f\T\i) the on-shell T-matrix. 

The results of Ref. [ID rely on the simple observation, 
that the on-shell matrix elements of the many-body T- 
matrix between single particle states, (pa\T\p' cr') , deter- 
mine both the total (<7t t) and the elastic (<J c i) scattering 
cross sections of the conduction electrons (or holes) at 
T = temperature. The total scattering cross section of 
an electron of momentum p and spin a is given by the 
optical theorem as 



°total 



Vf 



-lm(pa\T\pa) 



(5) 



where the velocity of the incoming electron is approxi- 
mated by the Fermi velocity, vp. In case of elastic scat- 
tering, an incoming single electron state is scattered into 



an outgoing single electron state, without inducing any 
spin or electron- hole excitation of the environment. The 
corresponding cross section can be expressed as 



1 
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dp' 
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with £ the energy of the electron measured from the 
Fermi surface. In contrast to <7°J, the total scatter- 
ing cross section also includes those scattering processes, 
where some excitations are left behind, and thus the out- 
going state is not a single particle state. The inelastic 
scattering cross section associated with these processes 
is thus clearly the difference of these two cross-sections: 



CT inol — CT total °el 



(7) 



These processes are schematically shown in Fig. [T] 

' .t 




FIG. 1: Sketch of (a) elastic and (b) inelastic scattering pro- 
cesses. In case of an inelastic scattering the outgoing electron 
leaves spin- and/or electron-hole excitations behind. 

For quantum impurities in a free electron gas it is more 
transparent to introduce angular momentum channels, 
L = (l,m), and define the scattering states in terms of 
radially propagating states |p, a) — > ||p|, L,a): 



Ip.o-) 



dn p Y L (p)\\p\,L,a) , 



(8) 



with Yl(p) the spherical functions, and g — f?(|p|) the 
density of states of the conduction electrons. In this basis 
the on-shell 5 and T-matrices become matrices in the 
angular momentum quantum numbers that depend only 
on the energy u> = £(p, a) of the incoming particle. The 
5-matrix can then be expanded as 

(pa\S\p'a') = -J^ Y l (P) 8 W^V M Yv (p') , (9) 

and the on-shell T-matrix is given by a similar expression, 
the coefficients of the expansions being related by 



SL.a- L'a> M = 5l,L>5<7,<j' + i th,a; V a' M 



(10) 



The eigenvalues 5a of the matrix s_l,o-; iv must all be 
within the complex unit circle for any to, and they are 
directly related to the inelastic scattering cross section. 
To see this, let us consider the case of scattering only 
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in the s-channel (L = 0), and assume spin conservation. 
In this case Sl !(T - iv becomes a simple number, s(u>) — 
1 + U(oj), 

s(ui) = 2ng (per|<S|p'o-') = 2ng S(ui) 

and the inelastic scattering cross section can be expressed 
as 

CinclM = -3- (1 - |s(^)| 2 ) 

= 4(2ImtH-|tH| 2 ), (11) 

where we assumed free electrons of dispersion £ = 
p 2 /2m — fi with a Fermi energy /i and a corresponding 
Fermi momentum pp. Eq. (jlip implies that the scatter- 
ing becomes totally elastic whenever s(u>) is on the unit 
circle, and it is maximally inelastic if the corresponding 
single particle matrix element of the S- matrix vanishes. 
The former situation occurs at Fermi liquid fixed points, 
while the latter case is realized, e.g., in the two-channel 
Kondo model or the two-impurity Kondo model. The to- 
tal scattering cross section, on the other hand, is related 
to the real part of s(u>) as 

2ir 2tt 
ototM = -0- (1 - Re{s(w)}) = -5- Im t(w) . (12) 
Pf Pf 

It is easy to generalize this result to the case of many scat- 
tering channels, and one finds that inelastic scattering 
can take place only if some of the eigenvalues of sl.o-- L'a' 
are not on the unit circled 

To compute the off-diagonal matrix element 
(per|7~|p'er') , we first relate it to the conduction 
electrons' Green function through the so-called reduc- 
tion formula detailed in Section HT1 28 

(p,a|T|pV)= (13) 

_ [ G °]±p,±a(0 G±p±<x,±p' ±a'(0 [G°] ±p , ±(J ,(£) . 

Here the ± signs correspond to electron and hole states 
with £ > and £ < and of energy E = |£|, Go denotes 
the free electron Green's function, and G the full many- 
body time-ordered electron Green's function. Thus the 
positive frequency part of the Green's function describes 
the scattering of electrons, while the negative frequency 
part that of holes. Strictly speaking, our derivation of 
this formula only holds for Fermi liquids, i.e., for models 
where the ground state has no internal degeneracy and 
can continuously be related to a non-interacting ground 
state. The consideration of internal ground state degen- 
eracy needs some care, and the definition of inelastic scat- 
tering in that case is not straightforward. 27 However, the 
finite temperature results of Ref. [24| show that Eqs. ([5]), 
©, and ([7]) together with (TIB")) provide the physically 
meaningful definition even in this case at T = 0. 

According to Eq. (fTB")) . to compute the inelastic and 
elastic scattering cross-sections, we need to evaluate the 
self-energy of the conduction electron's Green function. 
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FIG. 2: Renormalization group flow of the eigenvalues of the 
single particle S-matrix for the single- and two-channel Kondo 
models. For the single-channel Kondo model the scattering 
becomes elastic (\S\ = 1) both in the limit of very large and 
very small energies, ui — > 00, and ui — > 0, respectively. For the 
2-channel Kondo model s{uj — > 0) = 0, implying that half of 
the scattering remains inelastic even for ui — > 0. 

We do this by relating the self-energy to some local 
correlation function, that we then compute either ana- 
lytically within some approximation, or numerically us- 
ing the powerful machinery of numerical renormalization 
group (NRG) i2£ This step depends on the specific impu- 
rity model at hand, and can be achieved through equation 
of motion methods, 30 diagrammatically, 13 * or through a 
straightforward path integral treatment, as we do it in 
Section IIIII The recently-formulated scattering Bethe 
ansatz approach can possibly provide a way to avoid this 
numerical computation, and determine the full energy- 
dependence of the 5-matrix analytically 3 ^. 

Let us close this Introduction by illustrating the power 
of this formalism through the simple examples of the 
single- and two-channel Kondo models, defined through 
the Hamiltonian: 

/ 

H = ^ t y t £p a pcr,a a po\a 
a— 1 p,cr 

J - f 

a—1 p,p' 

an' 

Here a^ paa creates a conduction electron with momentum 
p, spin a in channel a, S = 1/2 is the impurity spin, and 
/ = 1 and / = 2 for the single- and two-channel Kondo 
models, respectively. 

In both models the T-matrix can be related to the 
Green's function of the so-called composite Fermion op- 
erator, F aa = p & ' v era 1 a-pcr' ar°- which can then be 
computed using NRG^S The evolution of the eigenvalue 
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of the numerically obtained S-matrix is shown in Fig. [2] 
In both cases, s(ui) = s(lo/Tk) is a universal function 
that depends only on the ratio uj/Tk, with Tk the so- 
called Kondo temperature, Tk ~ Ep e _1 / eJ , with Ep 
the Fermi energy and g the density of states at the Fermi 
energy for one spin direction . 32 ' 33 

For the single-channel Kondo model the scattering be- 
comes elastic both in the limit of very large and very 
small energies, u> ^> Tk, and u> <C Tk, respectively, 
where the eigenvalues lie on the unit circle. The rea- 
sons are different: At large energies conduction electrons 
do not interact with the impurity spin efficiently. At 
very small energies, on the other hand, the impurity's 
spin is screened and disappears from the problem apart 
from a residual phase shift of n/2 and an irrelevant local 
electron-electron interaction^ 4 - The maximum inelastic 
scattering is reached when the eigenvalue s(lo) is closest 
to the origin, i.e., at energies in the range of the Kondo 
temperature, u) rs Tk- 

For the 2-channel Kondo model, on the other hand, 
s(ui — > 0) = 0, implying that the scattering is maximally 
inelastic even at the Fermi energy, uj — > 0. This prop- 
erty of the S-matrix has been first noticed by Maldacena 
and Ludwig 3 ^, and is characteristic of a non-Fermi liquid, 
where incoming electrons do not scatter into an outgo- 
ing single electron state, even at the Fermi energy. Note, 
however, that the vanishing of the 5-matrix does not im- 
ply that all the scattering is fully inelastic. In fact, from 
Eqs. © and © it follows that 

_.2CK 
2CK _ „2CK _ "tot 
°incl — "el — 2 ' 

i.e., half of the scattering remains still elastic. In other 
words, in the elastic channel, the unscattered and scat- 
tered s-electron wave functions are completely out of 
phase, and therefore there is no net outgoing particle 
in the s-channel. 

The structure of the flow of s(u>) appears more directly 
in the in the energy dependence of the various scattering 
cross-sections shown in Fig. [3] for these two cases. In the 
single-channel case, it is quite remarkable that the low- 
energy ~ to 1 inelastic scattering cross section expected 
from Fermi liquid considerations 34 is limited only to the 
regime u> < 0.05 T K , and for 0.05 T K < u) < 0.5 T K 
the inelastic scattering cross section is quasi-linear. Fur- 
thermore, above Tk a wide plateau appears (rather then 
a peak), where (t(uj) is large and almost independent of 
the energy ui of the incoming particle. Both features 
appear also in a finite temperature calculation^ are in 
quantitative agreement with the experimental results of 
Refs. [j.22 on magnetically doped wires, and provide a 
possible explanation of the observed saturation of the 
dephasing time in some experiments on dephasing^. As 
we show in Section lVTl these universal features are robust 
and present in the Anderson model too. 

It is important to emphasize here that in the present 
paper we computed the inelastic scattering rate of elec- 
trons, rather than that of quasiparticles. This is moti- 
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FIG. 3: Top: Inelastic, elastic, and total scattering rates for 
the single-channel Kondo model in units of <to = An/pp for 
a spin up electron. o"i nc i scales approximately linearly with u> 
for 0.05 T K < (J < 0.5 Tk, behaves as ~ lo 2 for lu < 0.05 T K , 
and remarkably, shows a plateau above Tk- Bottom: Same 
for the two-channel Kondo problem for a spin up electron in 
channel 1. Note that the inelastic scattering cross section 
remains finite at ui — 0, and is precisely half of the elastic 
scattering cross section. 



vated by the trivial observation that in a real experiment, 
the external electromagnetic field couples with a minimal 
coupling to the bare conduction electrons. Precisely for 
this reason, the Kubo formula contains the conduction 
electron current operators, and also, the relevant quan- 
tity to determine dephasing is thus the inelastic scatter- 
ing rate of electrons. This is what we have computed 
here and that has been computed in Refsi^ 3 - and 2 ^. 

The definition of quasiparticles depends on the context 
in which they emerge. If one defines them as stable ele- 
mentary excitations of the vacuum, as Nozieres did 3 ^, or 
as they appear in Bethe ansata 4 ^ , then, by definition, 
these quasiparticles do not decay at all at T = and 
scatter only elastically 3 ^. 

Such quasiparticles are, however, usually complicated 
objects in terms of conduction electrons. For this reason 
they are typically not minimally coupled to the gauge 
field, and therefore, the current operator in the Kubo 
formula is a very complicated many-body vertex in the 
language of quasiparticles. Excepting for u> = 0, a real 
conduction electron is composed of many such stable 
quasiparticles, and it decays inelastically even at T = 
temperature, even if quasiparticles do not. In the Kondo 
model, at the Fermi energy quasiparticle states are simple 
phase shifted conduction electron states. However, the 
connection between quasiparticles and conduction elec- 
trons is not trivial for any finite energy. Therefore, if one 
considers inelastic scattering at a finite energy, one must 
precisely specify how finite energy quasiparticle states are 
defined, how they couple to a gauge field, and how a fi- 
nite energy electronic state is decomposed in terms of 
these quasiparticles. Unfortunately, except for the Bethe 
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ansatz, we are not aware of any work which would provide 
this necessary connection in sufficient detail, and would 
go beyond a simple heuristic treatment (which might 
still give the correct result). In the present framework, 
we avoided this difficulty by formulating the problem in 
terms of electrons rather than quasiparticles. 

The paper is organized as follows: In Section [IT] we 
present the derivation of the reduction formulas. In Sec- 
tion [TTT] we determine the T-matrix for the Kondo model. 
In Sections llVllVl and lVII wc present results on the inelas- 
tic scattering rate for the single- and two-channel Kondo 
model and for the Anderson model, respectively. In Sec- 
tion I VIII the results are summarized. In the Appendix 
some details of the derivation of the T-matrix for the 
Anderson model is discussed. 



II. REDUCTION FORMULAS 

A. Definition of scattering states in the Heisenberg 
picture 

Although reduction formulas are often used in the liter- 
ature in a heuristic way, apart from the derivation of Lan- 
greth for the Anderson model^ we do not know of any 
work that would establish a rigorous connection between 
the single-particle matrix elements of the T-matrix and 
the conduction electron's Green's function for a general 
quantum impurity problem. Here we therefore present a 
short derivation of the reduction formulas by generalizing 
the procedure used in the domain of particle physics i 28 ' 37 

In this section, following the field theoretical language, 
we shall use the Heisenberg picture, and describe scat- 
tering in terms of the field operators^- ip a (x), where 
we introduced the four- vector notation, x = (t, x). The 
evolution of this field operator is described by the time- 
dependent Hamiltonian, with the interactions switched 
on and off adiabatically with a rate 8 — > 0, 

H = H(t)=H Q + e- s ^H int . (15) 

Here Hq denotes the non-interacting Hamiltonian, 

H = J d 3 x - /#a(z), (16) 

with A the Laplace operator, and \x the chemical poten- 
tial of the electrons. The interaction part H ln t does not 
need to be specified at this point, and depends on the 
particular model considered. For the sake of simplicity 
we assume that the quantum impurity interacts with free 
electrons, but the procedure described can be generalized 
for electrons with more complicated dispersions, too. 

Within the Heisenberg picture, states are independent 
of time, and all non-trivial scattering is incorporated in 
the time evolution of the fields. Scattering states can be 
defined through the asymptotic form of the field oper- 
ators. Incoming and outgoing scattering states can be 
defined based on the simple observation that for times 



t — > ±oo the equation of motion of ipa(x) is generated by 
Hq, and therefore ip a {x) behaves asymptotically as a free 
field: 

tp a (x,t -oo) -> J -^^e- ip - x a pcr , in , (17) 

where a pCT ,i n = a po% i n are just the annihilation operators 
of incoming (one particle) scattering states. Here for the 
sake of compactness, we introduced the four-momentum, 
p = (£, p), with £ = £(p) = p 2 /2m — \i the energy of the 
conduction electrons, measured from the Fermi energy, 
and j? • x = £ i — px. The operators a pa ;m = ap<r,in 
satisfy standard anti-commutation relations: 

K*^,J = (2t) 3 <W £(P-P') • (18) 

Note that the operators a 1 a in do not create free elec- 
trons, rather, they are creation operators of incoming 
electrons in scattering states, which are asymptotically 
free^- 

The operators a p a in can be used to construct incoming 
single particle scattering states, |p, cr); n . For electrons, 
i.e., excitations of momenta larger than the Fermi mo- 
mentum, |p| > pp, these scattering states can be simply 
defined as^ 

|p, a) in = a+ ;CTiln |0) = lim / d 3 x e"** ^(x)|0) . 

(19) 

Outgoing single electron scattering states, \p, a) out = 
|p, c)out, can be defined in a similar way, by expanding 
the field ip a {x, t — > +oo), 

|p, a) out ee 4 >out |0) = lim / d 3 x e-*-« ^(x)\0) . 

(20) 

Incoming and outgoing hole states must be defined 
slightly differently, because an incoming hole of energy 
E > 0, momentum p, and spin a is created by removing 
an electron of energy £ = —E, momentum — p, and spin 
—a from the Fermi surface. In other words, incoming 
hole scattering states are defined for |p| < pp as 

|Pj f)in/out = a-p : -<r,in/out|0) (21) 

= lim / d 3 x e"^ iMaOlO) (|p| < p F ) . 

B. Reduction formulas and Green's functions 

We proceed to derive a general relation involving 
Green's functions to express the off-diagonal matrix ele- 
ments 

* <p ff|f|pV) = out (p «r|pV)ia (p?V) (22) 

for electronic excitations with |p| > pp first. Using the 
asymptotic expression Eq. (|19[) this matrix element can 
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be expressed as 



out(p c|p' a') ir 



lim „t(p,cr| / dye 
Integrating by part we obtain 



-«p'y 7 /,t 



V^,(y)|0) . (23) 



lim 



d 3 y e-^l(y) 



dyo 



PAv) 



a 



t 

p' : cr' ,out " 



The last term does not give a contribution to the matrix 
element for p ^ p' , therefore we drop it. The rest can be 
expressed as 



4 ° 



oyo 



VAv) 



dyo 



tfAv), 



where we obtained the r.h.s. of this equation by using 
the fact that p 1 is on the energy shell, and therefore 
p' e~ lp v — (— ^;A y — /i)e~ jp y , and then by integrating 
by part with respect to y. Thus the off-diagonal matrix 
elements of the ^-matrix simplify to 



out(P cr|p' a') in = 



i I d 4 ye~ v y 



(24) 

out<P^(y)|0) 



Using the asymptotic relation of the outgoing states, 
20| . we can now write the full matrix element as 



out(P trip' <?') 



i d ye 



-ip y 



. oyo 



x lim / d 3 xe^(0|^(x)^(?/)|0). 



(25) 



Once again, we convert the last integral into an integral 
over the whole space-time, which yields 



lim 

Xq — >00 



d 3 x e^<0|^(z)V^(j/)|0) 
d 



d'x—e^iOlT^ix^liym (26) 



where the time ordering operator T has been inserted 
to assure that the xq — + — oo contribution vanishes by 
Eq. (H|): \hn^^- oo fd 3 xe i **{0\1>l(y)il> tr (x)\0)=0. 

We can manipulate the remaining expression in the 
same way as before to finally obtain 



out(p <r\p' <r')iu = - I d*x d*y e^e** 



(-% ~d X0 - i?o(x)) (T^(x)^(y)) (i *8 yo - ffo(y))] , (27) 



where the arrows indicate forward and backward dif- 
ferentiation, respectively. Observing that the operator 
{—i d Xa — Hq(x.)) 5(x — x') is simply the matrix clement 
of the inverse of the non-interacting Green's function, 

{x'\G^\x) = S(x'-x) (-i^ Xo -H (x)) , 

= (i t x , o - ffo(x')) S(x' - x) , 

Eq. lj2"T|) can be simply expressed as 



out(P cr|p' a') 



(28) 



i / d 4 x d 4 y e ipx ~ ip y 



Gq * G * Gq 



x.y 



with '*' the four-dimensional convolution operator, and 
G the usual interacting Green's function, 



Gaa'{x,y) = -i {Tip^x^l^y)) . 



(29) 



The Fourier transformation of (|2"8|) then yields 

out{P cr\p' cr')in = (30) 

- i G ijp)G^,(p,p>)G ^(p>). (|p| > PF ) 

Translational invariance in time further implies 

G^fajJ) = 27r^(p')-C(p))G pCTipV ,(^(p)). 

Inserting this into (|30[) and comparing it with Eq. ^ 
yields Eq. (p)) for £ > 0. 

The derivation for holes follows exactly the same lines 
excepting that the matrix element to be computed is now 

out (p a\p' a') in = (0\at p ^ out a^ p> ^, . n \0}, (31) 

and correspondingly, the final expression of the 5-matrix 
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element now reads for p ^ p' 

out<P Cr\p' (j')in = 



(IpI <pf) 

-*(-p) ■ 



III. T-MATRIX OF THE KONDO MODEL 

For practical calculations, one needs to determine the 
T-matrix of Eq. (fT3")) somehow. For most really inter- 
esting cases this can be done analytically only approx- 
imately and in a limited energy range, and numerical 
methods must be used. The most adequate way to per- 
form the calculation is to first relate the T-matrix to some 
local correlation function that can then be computed us- 
ing Wilson's numerical renormalization group method. 29 
To establish the desired relation, one can use equation of 
motion methods^ or do diagrammatic perturbation the- 
ory and sum up the diagrams up to infinite order^ but 
here we show yet another rather elegant way, in terms of 
path integrals. 

Although this method works for essentially any quan- 
tum impurity problem, here we show how it works for the 
Kondo model, already defined in the introduction (see 
Eq. HH)). The application of this method to the Ander- 
son model is discussed in the Appendix. To use a field 
theoretical formalism, following Abrikosov, we represent 
the impurity spin S by fermionic operators, 



^ ^ f erf a 



1 



(32) 



Next, we define the generating functional for the conduc- 
tion electron Green's functions as follows: 



Z[ Va ,r] a ] = / D [a a a a ] D [f a f a ] e~ l s e 



i S irf-a+ia-rj 



where the tilde on the second integration measure indi- 
cates that one must impose the constraint f\f a = 1 
when performing the path integral, and we introduced 
the shorthand notation: rj - a = a J dtrj p a (t)a PtCr (t) . 
The action S in Eq. (f3"3"|) consists of three terms, S = 
S e + Sf + Si n t: The first term, S e describes the conduc- 
tion electrons, 



S e = -J2 f dtdt'a pa {t)[G%\t-t')a pa {t') (34) 

p,a J 



with G° the time-ordered free electron Green's function. 
The term Sf = —i J2 a J dt f a (t) 4tfa(t) generates the 
spin dynamics, while the last term simply describes the 
interaction: 



Shit = — ^2 J dt ^(*) <Vr(*)°W a pV'(*) • ( 35 ) 
p-p' 

<t,<t' 

The full time-ordered Green's function is related to Z by 



G 



po\p a 



(*-*') = 



1 <5 2 lnZ [77,77] 



(36) 



rj— rj— 



We can derive the required identity by simply shifting 
the integration variable in Eq. (I33|) , 



i P *(t) -> a prT (t) - J dt' G° p (t - t') r) pa {t') . (37) 



(33) As a result, the exponent in Eq. ([33]) transforms to: 



S-rj-a-a-ri S + £ / dt dt' rj p „(t)G p(r (t - t')i lpa {t') - dt dt' (F a (t)G° p<7 (t - tf) Vp<7 (t') + h.c.) 

+ Jj2jdt dt' dt" v pa (t)G° p At - t')S{t')a aa ,G% a {t' - t") Vp , a (t") , 
p,p' 

CT.C' 

where we introduced the composite fermion field, F a (t) = Y] p ^ S(t)a acT >a p . a >(t). Carrying now out the functional 
derivation of (|36p we obtain the following simple relation 

G pa>p , a , (t-f) = S p , p , G° pcr (t-t') + J dtdt' G° pa (t-t) Ui-F) l(S)a a „, - i^{F„(t)F G° p , a ,(t' -t'). (38) 

I 

The average in this expression must be carried out by the corresponding time-ordered Green's function. Corn- 
computing the appropriate path integral, and results in paring Eqs. (|3"5|) and (p~5|) . and using the analytical prop- 
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erties of the time-ordered and retarded Green's functions at T — temperature, we finally obtain the relations: 



ln\{T Ta . T<T ,{uj)} = 7T— Q F , aa '{ijS) , 

-rRe{T T ^H} = ~(S)<w + ^ 4 ckj' gF ^' ( ^ 
z 4 / — 



(39) 



T = 



+ 



+ 



Q 



+ X^A- 



+ 



T 



FIG. 4: Diagrammatic representation of the T-matrix in the 
Kondo problem. Dashed lines denote pseudo-fermion propa- 
gators and describe the evolution of the impurity spin, while 
continuous lines denote free conduction electron propagators. 
Filled circles stand for the exchange interaction J. The first 
term of the T-matrix is simply proportional to the expecta- 
tion value of the impurity spin, it is frequency independent, 
and vanishes in the absence of magnetic field. The second 
term can be identified as the composite fermions correlation 
function. 



with QF,aa'{u) the spectral function of the composite 
Fermion's Green function, and r = sgn £. Note that 
scattering takes place only in the s-channel, and there- 
fore these matrix elements do not depend on the incoming 
and outgoing momenta of the excitations. Eqs. (|38[) and 
(|39[) can be easily visualized in terms of diagrammatic 
perturbation theory, as shown in Fig. |4j The spectral 
function appearing in Eq. (|39p is just a local correlation 
function that can be easily obtained through the NRG 
method. 29 In the following sections we shall primarily use 
this method to compute the single particle T-matrix and 
the inelastic scattering rates of the basic quantum im- 
purity models, the single-channel Kondo model, the two- 
channel Kondo model, and the Anderson model. Calcu- 
lations for the spin S = 1 Anderson model have been 
performed in Ref. [H 



IV. INELASTIC SCATTERING IN THE KONDO 
MODEL 

In the previous sections we have related the single par- 
ticle T-matrix and therefore the elastic and inelastic scat- 
tering amplitude of electrons with local correlation func- 
tions through the reduction formulas. In this section we 
shall use these results to analyze the T = temperature 
scattering properties of the Kondo model using Wilson's 
NRG 29 . However, before presenting detailed numerical 
results, let us shortly discuss what one can learn from 
simple perturbation theory. 

Let us discuss the high-energy scattering of conduction 
electrons in the absence of external magnetic field. In this 
limit one can attempt to do perturbation theory, and in 
first non-vanishing order one obtains 



t(u>) — i 2kq T(u>) 



-i— S(S + 1) j + 



(40) 



where the dimensionless coupling j = gj has been in- 
troduced. Summing up the leading logarithmic diagrams 
amounts in replacing J by the renormalized coupling, and 
gives 

t{w » T K ) » *^S(S + 1) 1 , 
I In (uo/Tk) 

with T K ~ E F e~ x l jQ the Kondo temperature. Thus, 
in leading logarithmic order, the total scattering cross 
section is given by: 



1) 



1 



Pf In {u/T K ) 

The first non-vanishing contribution to the elastic scat- 
tering cross section, on the other hand, scales as o~ c \ ~ 
\t\ 2 ~ j 4 , and therefore a e i(u) asymptotically behaves as 



<7 cl (a; » T*:) « _ S 2 (S+l)\ , 



hyV/Tx) 



This implies that asymptotically, all the scattering is in- 
elastic 



d(w » T K ) « <7totM w -r S ( S + 1) 



1 



p% \n\u/T K ) 



(41) 
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This is a very surprising result, and contradicts to the 
conventional wisdom, which tries to associate inelastic 
scattering with spin-flip scattering from a free spin. In 
fact, this rather non-trivial result has been explained in 
Ref. [3^ in the following way: At high energies, incoming 
electrons are scattered by the impurity spin fluctuations. 
These fluctuations can absorb an energy of the order of 
~ Tk, and therefore the energy of the incoming electron 
is not conserved in leading order, but it typically changes 
by a tiny amount, Slo ~ Tk- In the most pedestrian per- 
turbative approach this tiny energy transfer is neglected 
and therefore one concludes incorrectly that the energy 
is conserved in leading order. 

We can also relate the cross sections above to scatter- 
ing rates. Assuming a finite concentration ni mp of mag- 
netic impurities, we can compute the impurity averaged 
conduction electron Green's functions and from that the 
conduction electron lifetime: 



tip VF CTtot(w) W 



tt5(5 + 1) 
2 g \n 2 {uj/T K ) 



(42) 



In fact, the first part of this equation gives a general rule 
to connect various cross sections to the corresponding 
scattering times, and for the inelastic scattering rate, e.g., 
we have 



1 

Tin el 



ni m p VF CTincl(^) • 



(43) 



For very large frequencies, again, the inelastic scattering 
rate is approximately equal to the elastic scattering rate: 



1 

Tin el 



ttS(S+ 1) 
2 g ln 2 {u/T K ) 



(44) 



Note that this rate is a factor of 3/2 larger than the 
Nagaoka-Suhl expression, which only takes into account 
spin flip processes^ 

For energies \lu\ <C Tk perturbation theory breaks 
down, and it is more appropriate to use Nozieres' Fermi 
liquid theory, which states that at the Fermi energy scat- 
tering is completely elastic, and^i 



t a {uj = 0+) = 2ng T^ ck {uj = 0+) 



2sin<L e lS ° 



(45) 



where we now allowed for a magnetic field pointing along 
the z-direction, and 5 a stands for the phase shifts of elec- 
trons with spin a at the Fermi energy. Eq. (|45|) then 
yields 



P'f 

Cmchcr^ — > 0) = 



4?r • its 
— sin (d 



(46) 
(47) 



The maximum total scattering cross section is reached in 
the unitary limit, 8 a = ir/2. 

Let us now proceed and compute the various scatter- 
ing cross section using Wilson's NRG—. In the previous 
section, we showed how the imaginary part of the single 
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FIG. 5: Renormalization group flow of the eigenvalues of the 
single particle .S'-matrix for the single-channel Kondo model 
in presence of a local magnetic field. 



particle T-matrix is related to composite Fermion's spec- 
tral function, Qp{oS). Within NRG, spectral functions of 
local operators are computed using their Lehman repre- 
sentation. The imaginary part of the T-matrix, related 
to the total scattering rate has already been computed 
in this way by Costi to obtain the magneto-resistivity of 
Kondo alloys^ To evaluate the inelastic scattering am- 
plitude, however, one needs to go one step further and 
compute the real part of the T-matrix as well through 
a Hilbert transformation, Eq. (f3"9"]l . In such a calcula- 
tion it is essential to have high quality data. The most 



challenging task is to obtain the correct 



low energy 



behavior of the inelastic amplitude since we get this small 
quantity as a difference of two quantities of the order of 
unity. Therefore it is also crucial to get the normaliza- 
tion factor of T correctly. In case of the single-channel 
Kondo problem this can be obtained through the Fermi 
liquid relation, (|45[) : This relation connects the normal- 
ization of t(u>) to the phase shifts at the Fermi energy, 
which we extract from the NRG finite size spectrum very 
accurately^. 

The renormalization group flow of the eigenvalues of 
the single particle S-matrix has already been shortly dis- 
cussed in the introduction in the absence of magnetic 
fields (see Figfl]). The eigenvalues s(ui) of the S'-matrix 
lie within the complex unit circle, and inelastic processes 
are allowed only when the eigenvalue s(w) of the S'-matrix 
satisfies |s(w)| < 1. Incoming particles of high enough 
energy (lj — > 00) do not see the impurity, and therefore 
s(lo) — > 1, corresponding to the weak coupling fixed point 
with phase shift 6 sw 0. As expected for a Fermi liquids, 
at the Fermi energy (lj — ► 0) the s(lj) approaches the unit 
circle again, s(u> = 0) = — 1, corresponding to the strong 
coupling fixed point of the Kondo model characterized 
by a phase shift 6 — tt/2. 
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Application of a local magnetic field makes the flow 
more complicated (see FigJSJ): At low energy the sys- 
tem still behaves like a Fermi liquid but the position of 
the point where the s(o>) approaches the unit circle now 
varies with magnetic field. This is due to the magnetic 
field dependence of the scattering phase shifts at zero 
frequency. 

For intermediate energy values of the incoming elec- 
tron \s(u>)\ < 1 and the inelastic scattering cross-section 
is non-zero. The total, elastic and inelastic scattering 
cross sections of an electron scattered off a magnetic im- 
purity are shown in Figj6] As expected, the inelastic am- 
plitude always vanishes at the Fermi level. In the lower 
panel we show the inelastic scattering rate as compared 
to the Nagaoka-Suhl formula. The numerical results are 
consistent with the analytical expression (|4Tj) at large 
energies, while for energies much smaller than Tk we re- 
cover the quadratically vanishing inelastic rate expected 
from Fermi liquid theory. 42 Note that the Nagaoka-Suhl 
approximation systematically underestimates the inelas- 
tic scattering rate by a factor of 2/3 since it considers 
any spin-diagonal process as elastic scattering. At high 
energies, however, in leading order all the scatterings are 
inelastic since even a spin diagonal process breaks up 
the Kondo singlet and leaves the system in an excited 
state, and therefore it cannot be elastic. Apart from this 
prefactor, the Nagaoka-Suhl result is perfect at high en- 
ergies, however, it starts to deviate strongly from the 
numerically exact curve at approximately IOTr-, and it 
completely fails below the Kondo temperature Tk- At 
energies well above Tk almost all the scattering is inelas- 
tic, i.e. the inelastic amplitude varies as ~ In - (ui/Tk) 
while the elastic part vanishes faster as ~ In - (uj/Tk), 
in agreement with the analytical results. 

Even though the numerics recover the expected asymp- 
totics, interesting features appear both in the low and 
high energy part of the scattering properties. First, as 
shown in Fig{3j the ai ne i ~ ui 2 regime appears only at en- 
ergies well below the Kondo temperature, and we find 
that the inelastic scattering rate is roughly linear be- 
tween 0.05Tk < u) < 0.5Tk- Even though our calcu- 
lation is done at T — temperature, we expect that 
&iae\(T, to = 0) behaves very similarly to a- mc i(T = Q,u>). 
Our results are thus consistent with the existing exper- 
imental data, explain the linear behavior observed in 
many experiments;^^ and surprisingly even quantita- 
tively fit the finite temperature experimental curves^ 
Of course, in reality a finite temperature calculation is 
needed which has been performed in Refill. 

Another remarkable feature is the broad plateau in the 
inelastic scattering cross section above the Kondo scale, 
where the energy-dependence of the inelastic scattering 
rate turns out to be extremely weak. This weak energy- 
dependence provides a natural explanation for the exper- 
imentally observed plateau of the dephasing rate in many 
experiments.— £ 

The inelastic scattering amplitude in presence of a 
magnetic field is shown in Fig(7] Applying a local mag- 
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FIG. 6: Upper panel: Inelastic, elastic, and total scattering 
rates for the single-channel Kondo model in units of ao = 
4-7r/p|n at T = and B — 0, as a function of the incoming 
electron's energy. Only the electronic contribution (u > 0) is 
plotted. For u) 3> Tk the inelastic rate as well as the total 
scattering rate varies as In -2 (Tk /lo) while the elastic part 
decays as ln -4 (TR-/u;). The lower panel shows the ~ to 2 and 
In -2 (Tk /w) regimes for w <C Tk and u> ^> Tk, respectively. 
The Nagaoka-Suhl approximation, corrected with a factor 3/2 
(see the text) is also shown. 



netic field breaks the spin symmetry of the scattering and 
changes the inelastic scattering properties of spin up and 
down particles dramatically. Already a relatively small 
magnetic field B ~ O.ITr- results in a very strong spin 
asymmetry of the inelastic scattering. For B ~ Tk the 
effect is even more dramatic: At this field the impurity 
is practically polarized and aligned with the direction of 
the field and points upwards. As a result, an incoming 
spin up particle cannot flip the local spin in a first order 
process, and higher order processes are needed to gener- 
ate inelastic scattering. A spin down electron or hole, on 
the other hand, can exchange its spin with the magnetic 
impurity, resulting in the maximum of the inelastic rate 
at energy w B and a very broad inelastic background for 

LO > B. 



V. INELASTIC SCATTERING IN THE 
TWO-CHANNEL KONDO MODEL 

In this section we shall present results for the two- 
channel Kondo model, the prototype of all non- Fermi liq- 
uid impurity mo dels j 43 1 44 In the channel symmetric case 
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FIG. 8: Renormalization group flow of the eigenvalues of the 
single particle S'-matrix for the two-channel Kondo model for 
various values of channel asymmetry. The curves correspond 
to ji +j 2 = 0.12. 



FIG. 7: Energy dependence of spin-dependent elastic and in- 
elastic scattering rates for the single-channel Kondo model in 
units of (To = A-K/pp, at T — and in the presence of a local 
magnetic field B. Note that the positive frequency side corre- 
sponds to spin down electrons, while the negative frequency 
side is for spin up holes. The impurity spin is polarized up- 
wards. 



there is two types of conduction electrons that try to 
screen the impurity independently leading to the over- 
screening of the local moment. This frustration of the 
screening processes manifests itself in the formation of 
a strongly correlated state which cannot be described in 
the framework of Fermi liquid theory. This unusual cor- 
related state manifests in the nonzero residual entropy, 
the logarithmic divergence of the impurity susceptibility 
and the power law behavior of transport properties with 
fractional exponents. 

Since the non-Fermi liquid behavior is a direct con- 
sequence of the frustration of the screening processes, 
any infinitesimal asymmetry in the couplings A = (ji — 
3i)l (ji +J2) leads to the appearance of another low tem- 
perature energy scale T* cx A 2 /Tk at which the sys- 
tem crosses over to a Fermi liquid: Electrons being more 
strongly coupled to the impurity form a usual Kondo 
singlet with the impurity spin, while the other electron 
channel becomes completely decoupled from the spin. 

In the 2-channel Kondo case, unfortunately, no Fermi- 
liquid relations similar to Eq. (|45|) are available. How- 
ever, there is an exact theorem by Maldacena and Lud- 
wig, that allows us to get the right normalization of the 
T-matrix. This theorem states that, at the two-channel 
Kondo fixed point, the single-particle elements of the S'- 
matrix vanish at the Fermi energy, S2Ck(^> — > 0) = 0^ 



and as a consequence 

hcniu = + ) = -i ■ 

This relation allows us to obtain the proper normaliza- 
tion of the numerically computed T-matrix even at the 
non- Fermi liquid fixed point. However, it also leads to 
the surprising result mentioned already in the introduc- 
tion, that exactly half of the scattering is inelastic at 
the Fermi energy, while the other half of it is inelastic. 
This counter-intuitive result can be understood as fol- 
lows: The identically zero single particle S'-matrix indi- 
cates that an incoming electron cannot be detected as 
one electron after the scattering event, and it "decays" 
into many electron-hole pairs. To get such a "decay", 
the scattering process must have an elastic component 
too which interferes destructively with the not scattered 
direct wave and results in the absence of the outgoing 
single particle amplitude in the s-channel. 

The universal flow of the eigenvalue of the S-matrix 
was shown in Fig. [2j In Fig|8]we show what happens if we 
make the couplings in the two channels slightly asymmet- 
ric. For any small asymmetry the Fermi liquid behavior 
reappears: The S-matrix in the more strongly coupled 
channel flows first close to the two-channel Kondo fixed 
point with s(oj) ~ 0, and then below the Fermi liquid 
scale ui < T* it suddenly crosses over to the strong cou- 
pling fixed point characterized with phase shifts 5 = n/2 
Similarly, s(uj) in the other channel also approaches the 
two-channel Kondo fixed point, but then it becomes sud- 
denly decoupled and therefore S flows to the s(u>) = 1 
fixed point. 

The inelastic scattering rates for the two-channel 
Kondo model are shown in Fig. [5] as a function of the en- 
ergy of the incoming particle. In the channel-symmetric 
case inelastic processes are allowed even at uj = 0, which 
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FIG. 9: Energy dependence of elastic and inelastic scattering 
rates for the two-channel Kondo model in channel "1" in units 
of ctq = 4-7T /p% , at T = and for different channel asymmetry, 
A = (Ji— JaJ/O'i+Ja)- The curves correspond to ji +J2 = 0.3. 



is a clear signature of the non-Fermi liquid behavior. The 
non-Fermi liquid nature is also reflected in the ~ y/uj sin- 
gularity of the scattering cross sections at u = 0. Note 
that this cusp is much less pronounced in the inelastic 
scattering rate. 

For A > the total scattering rate approaches the 
unitary limit in channel "1" below the Fermi liquid scale 
T*. For A < 0, on the other hand, the total scattering 
rate goes to in channel "1" below the Fermi liquid scale 
T*. In both cases, the inelastic scattering freezes out, and 
<7inei(w) shows a dip below T*, and it ultimately scales 
to as 



CTinclH « CSt. —5 • 



Note that the inelastic scattering cross-section is very 
similar for A > and A < 0, while the total scattering 
contributions differs dramatically in these two cases. 



VI. INELASTIC SCATTERING IN THE 
ANDERSON MODEL 

As a final example, let us consider the Anderson model 
defined by the Hamiltonian, 

pa a 

+ Ud\d ] d\d l (4<A + h.c.) , (48) 

(T,p 

where now d a denotes the local d-level's annihilation op- 
erator, U is the on-site Coulomb repulsion, and the con- 
duction band and the local electronic level are hybridized 
byV. 

The T-matrix for the Anderson model can be related to 
the d-level's Green's function, as first discussed by Lan- 
grethi^ The required relation can be trivially established 
the the path integral formalism presented in Section [Ml 
The final result of this derivation, which is to some extent 
discussed in Appendix lAl can be written as: 



Im{7^(w)} = irV g d , T a{^) 



V 2 / dJ 



(49) 



with Qd,aiS) the spectral function of the d-Fermion's 
spectral function, and r = sgn £. 

The ground state of the Anderson model is of a Fermi 
liquid. Therefore, the Fermi liquid relations (|45]) can be 
used again to properly normalize the T-matrix. As we 
discussed in Ref. [l3|, the Fermi liquid relations also imply 
that at the Fermi energy the eigenvalue of the single par- 
ticle iS-matrix lies on the unit circle, and therefore the 
inelastic scattering rate vanishes. 

The flow of the eigenvalue s(u>) is shown in Fig. [TO] 
This flow diagram is very similar to that of the Kondo 
model at low energies, however, a new interesting feature 
appears at to ~ U, where s(u>) displays a hook. This 
hook corresponds to largely inelastic scattering processes, 
which are associated with the charge fluctuations of the d- 
level. It is adequate to mention here that the low energy 
flow is not completely identical to that shown in Fig [5] 
for the Kondo model. The reason is purely technical: 
For the Anderson model we were using the self-energy 
trick invented by Bulla and coworkers^ to obtain higher 
quality results, and computed the u> ~ Tk part of the 
T-matrix (especially its real part) much more accurately. 

These features also appear in the various scattering 
rates, shown in Fig. [TT] for the symmetrical Anderson 
model with = —U/2. There we show the inelastic 
scattering rate for various ratios of U/A, A = ttqV 2 be- 
ing the width of the resonance. For moderate values of 
U I A the effects of U are minor in the total and the elas- 
tic scattering rate, however, rather surprisingly, one can 
see a clear maximum in the inelastic scattering rate at 
energies u> ~ U even in this case. Increasing U/A, the 
Fermi liquid regime and the charging regime separate, 
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FIG. 10: Flow of the eigenvalue of the S'-matrix for the 
Anderson model for A = 0.035, U = 0.04, 0.08, 0.16, 0.26, 
0.32, 0.40, 0.64 and e d = -17/2. 
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FIG. 11: Elastic, inelastic and total scattering cross section 
for the symmetric Anderson model £d = —(7/2, for various 
ratios of U/A. For moderate values of U/A the effects of U 
are minor in the total and the elastic scattering rate, while 
for large values of U / A the various scattering rates for ui <C U 
follow very nicely the behavior found for the Kondo model 



and two distinct peaks appear, now even in the total 
and elastic scattering rates. For large values of U/A 
the various scattering rates follow very nicely the behav- 
ior found for the Kondo model at low energies, and for 
Tk <C co <C U the elastic and inelastic contributions scale 
as ~ l/ln 4 (o;/r^) and ~ 1 / In 2 (uj/Tjc), respectively. It 
is remarkable, that the Hubbard peak at u> ~ U is essen- 
tially entirely inelastic. 

Fig. [T2] shows the same behavior on a linear scale for 




FIG. 12: Elastic, inelastic and total scattering cross section 
for the asymmetric Anderson model. The low energy part of 
the curves is again very similar to the one obtained for the 
Kondo model. 



the asymmetrical Anderson model with moderate inter- 
action strength. The low-energy part of the figure is 
again strikingly similar to the one obtained for the Kondo 
model. This is not very surprising, since the Kondo 
model is just the effective model of the Anderson model 
in the limit of large U/A and u> <C U, where charge fluc- 
tuations occur only virtually. It is remarkable that the 
quasi-linear behavior of t7i no i and the plateau are already 
present for these moderate values of U/A. 



VII. CONCLUSION 

In this paper, we discussed in detail the theory of in- 
elastic scattering from quantum impurities at T = tem- 
perature, as formulated in Ref. [T3 . and applied this for- 
malism to various cases. We computed numerically the 
flow of the S'-matrix eigenvalues s(u>) for three prototyp- 
ical examples of quantum impurity models, the Kondo 
model, the two-channel Kondo model, and the Ander- 
son model. As we discussed, inelastic scattering appears, 
once \s(u>)\ < 1, and the crucial difference between Fermi 
liquid models and non-Fermi liquid models is that for 
non- Fermi liquid models |s(o>)| < 1 even at the Fermi 
energy, to — > 0, while for Fermi liquids \s(to = 0)| = 1. 

We also determined the inelastic scattering cross sec- 
tion, (Ti n ei(w), for all these models. For the Kondo 
model and the Anderson model in the Kondo regime, 
i.e., for large interaction values, the low-energy part of 
Cinei(w) has features essentially identical to those of the 
Kondo model: Deep in the Fermi liquid regime one has 
cr in eiH - (oj/T k ) 2 , while for 0.05T K < u < 0.5T K a 
quasi-linear regime appears, above which tTi ne i exhibits a 
plateau with tJi ne i(a;) « est over a wide frequency range. 
These features are quite robust, and survive even in the 
case of electron-hole symmetry breaking. 
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We also find that at large frequencies the scattering 
becomes asymptotically inelastic, and the inelastic scat- 
tering rate scales as 



1 

Tin el 



tt5(5+ 1) 
2 g ln 2 (u/T K ) 



while the elastic scattering rate falls off much more 
rapidly as 



1 

Tel 



''imp 



tt 3 5(5 + 1) 
8 q ln 4 (io/T K ) 



This result implies that -contrary to common wisdom- 
even spin-diagonal scattering is inelastic at high ener- 

gieSiai 

In addition to these remarkable low-energy features, 
the Anderson Hamiltonian exhibits another, very inter- 
esting inelastic scattering peak at ui ~ U that corre- 
sponds to charge excitations. Rather surprisingly, this 
peak is present even in the weak coupling regime, where 
no Hubbard peak can be seen in the total scattering cross 
section. In the Kondo regime, on the other hand, this 
peak is essentially identical to the Hubbard peak that 
appears in the total scattering cross-section, and which 
corresponds to almost completely inelastic scattering. 

In the two-channel Kondo model, the prototype of all 
non-Fermi liquid models, inelastic scattering remains fi- 
nite even if uj — * 0, and is exactly half of the total scatter- 
ing rate. However, the tiniest channel symmetry breaking 
destroys this non-Fermi liquid state, and generates a new 
Fermi liquid scale, T*, below which inelastic scattering 
freezes out, and the scattering becomes totally elastic. 

The inelastic scattering rates computed here for the 
Kondo and Anderson models and their finite temperature 
versions computed in Ref. [24| are in quantitative agree- 
ment with recent experimental studies on magnetically 
doped mesoscopic wires excepting the limit of very small 
temperature, where a small residual inelastic scattering 
rate seems to be present . 21 i 22 i 23 The origin of this small 
residual inelastic scattering rate is not clear yet, it might 
be due to some structural defects caused by the implan- 
tation process, or just some magnetic ions located at the 
interface of the wire. The agreement is even more surpris- 
ing, since in reality, magnetic impurities are not of spin 
5 = 1/2 character, but have a rather complicated ci-level 
structure*^. They thus usually have a large spin associ- 
ated with them (typically 5 = 2 or 5 = 5/2 for for Fe, 
Cr, or Mn) subject to crystal fields, that does not couple 
through a simple exchange interaction to the conduction 
electrons. In reality scattering thus takes place in some 
d-electron channels. For 5 = 5/2, e.g., the Fermi liquid 
state forms due to screening in five d-channels. Unfor- 
tunately, these realistic impurity models are out of reach 
for NRG computations. 

In case of d- impurities, scattering cross-sections be- 
come also larger due to the many angular momentum 
channels that are open to scattering. Assuming spheri- 
cal symmetry, e.g., the angle averaged total and elastic 



scattering cross sections become 

ftot = -?y^ Im{t L (u)} 

P F r. 



(50) 
(51) 



i.e., the total, elastic and inelastic scattering cross- 
sections are about five times larger for ci-wave scatter- 
ing than for s-wave scattering considered in the usual 
Kondo problem. This must also be taken into account 
when computing the amplitude of the observed Kondo 
anomaly or that of the inelastic scattering rate. Finally, 
band structure effects may also play an important role 
in real materials, where the Fermi surface is not spher- 
ical, and the Fermi velocity depends on the direction of 
incidence 25 . 
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APPENDIX A: FIELD-THEORETIC 
DERIVATION OF THE T-MATRIX FOR THE 
ANDERSON MODEL 

Here we derive the T-matrix for the Anderson model 
following the lines of Sec. IIIII We first introduce the 
generating functional for the Green's functions 



Z(r),rj) = / D[a tT a a ]D[d <T d <T ]e- iS e i7ia+iria , (Al) 

where the source terms are defined as in Sec. IIIII The 
action 5 = 5 e + 5<j + Sh y b + 5j„ t consists of four distinct 
parts, with 5 e defined by Eq. (|3~4")) is identical to the one 
given in Sec lIIII and the remaining parts given by 



Sd = I dtdt'd a (t)[G Q d ] a {t-t')d a {t') , 

Shyb = / dt(a a {t,0)d <T (t) +c.c) , 

Sint = U I dtd^{t)di{t)di(t)d^{t) , (A2) 



respectively. Shifting the integration variables generates 
the terms 
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- rj-a-va^J2 I dtdt'rj p>a G° p(T (t - t') Vp , a - f dtdt'G° pa (t - t 1 ) (%Jt')d a (t) + c.c.) , 

(A3) 

and which, after functional differentiation with respect to ?7 P ,er, give rise to the identity 

<W«r/(t - = <5p jP /^,a'^(t - - iV 2 ^/ / dtdi'G° pa (t - t)(d«(t)d„,(t'))G p , a ,(t' - t'). 



(A4) 

The Fourier transform of this equation yields Eq. (|49| . 
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